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In this paper, we present a framework for getting a series of exact vacuum 
solutions to the Einstein equation. This procedure of resolution is based on a 
canonical form of the metric. According to this procedure, the Einstein equation 
can be reduced to some 2-dimensional Laplace-like equations or rotation and 
divergence equations, which are much convenient for the resolution. 
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1 Introduction 



is a highly nonlinear system of partial differential equations satisfied by the metric g^, 
where R^ v stands for the Ricci tensor and R the scalar curvature of the spacetime 
manifold defined respectively by 
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in which is the Christoffel symbol given by 
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is the energy-momentum tensor of matter, and the Greek indexes go from to 3 
or from t, z, x to y. In the vacuum domain we have 



T, 
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To get the exact vacuum solution of the Einstein equation is an arduous work[l, 2]. 
The conventional method of resolution is to analyze the symmetry of the metric, i.e. 
analyze the Killing vectors of the spacetime. Almost all the well-known solutions 
such as the Friedmann- Robert son- Walker metric, Bianchi universe, Lemaitre-de Sitter 
universe, as well as Schwarzschild metric and Kerr metric, Taub-NUT solution[2, 3, 4], 
are based on the symmetry of spacetime. 

Obviously, a good choice for the coordinate system of the spacetime manifold should 
be very helpful for analyzing the Einstein equation. The traditional choices are the 
Gaussian coordinate system or the harmonic coordinate system [5]. The former gives 
locally the following metric 
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while, the latter satisfies the following coordinate condition 
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However (1.6) and (1.7) are more convenient only for some theoretical analysis rather 
than for practical resolution. 

Fortunately, the problem can be much simplified in a specially chosen coordinate 
system. Under some natural assumptions we have found that the metric can be trans- 
formed into the following canonical form 
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with 

det(^) = -v 2 ab (1.9) 

and the inverse 
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where u,v,p,q,a,b are smooth functions of the coordinates (t,z,x,y) (See [6] for the 
proof of this canonical form). 

The canonical form (1.8) is an Alibaba's conjuration for solving the Einstein equa- 
tion. In what follows, we will show how to use this metric to solve the vacuum Einstein 
equation. 



2 Procedure of Resolution 



For the canonical metric (1.8), computing the Einstein tensor G^, we get 
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where V = V(t, x, y) is a function to be determined. 

In order that the following recursive procedure goes on, for the exact solutions, we 
must adopt an ansatz for a and b to make (2.2) integrable. Different ansatz leads to 
different family of solutions. In what follows we take the separating variable ansatz as 
a representative example to show how the procedure works, that is, we assume 



a = L(t, x, y)K(t, z) 



b = N(t,x,y)K(t,z) 



(2.3) 



where d z K ^ 0. 

Substituting (2.3) into (2.2), it is easy to see that 
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Using (2.3) and (2.4), we have 
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Solving G zx = with respect to p, we get 



where P, A are functions of t, x, y to be determined. 
Similarly, we have 

1 ( d 2 z q d 2 z Kd z q d z K n r d z K . 
Solving = with respect to q, we get 



where Q, B are functions of t, x, y to be determined. 

Substituting (2.3), (2.4), (2.6) and (2.8) into G xy = 0, we get 

d x Q - d y P = -P§£ + Q^f (2.9) 

and 

A = 5 = 0. (2.10) 

By G xx = or G yy = 0, we obtain 
d x P d y Q P (8 X L d x N\ Q (d y L d y N\ 1 (d t L d t N" 
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Therefore, taking (P, Q) as a vector in subspace (x,y), (2.9) and (2.11) provide the 
corresponding covariant rotation and divergence equations respectively. 
By G tz = 0, we get 

u = ~ (t + af) a ' 2 + 2a ' (yK) + WK + ^ + vr> (2 - 12) 

where a is a function independent of z to be determined. U (t, x, y) and W(t, x, y) 
satisfy 

= ^ 1 -gy (2.13) 

with 
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in which the covariant Laplace and gradient operators of a scalar in the subspace (x, y) 
are defined respectively by 
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By G tx = G ty = 0, we find a = a(t). We should point out that, G tx = G ty = also 
provide two linear Laplace-like equations of P and Q, which are consequences of (2.9) 
and (2.11). 

For G u = 0, a direct calculation leads to 
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Thus, we can take U, V, W, P, Q, a as unknown functions satisfying the almost 
linear partial differential equations (2.18), (2.13), (2.19), (2.9), (2.11) and (2.15), while, 
take K(t,z), L(t,x,y), N(t,x,y) as functions determined by coordinates, boundary 
conditions and initial conditions. Any one of the equations (2.18), (2.13) and (2.19) is 
quite similar to the Ernst equation [7, 8]. 

Furthermore, we can simplify the metric by setting r = K as new coordinate to 
replace z as follows. The line element of the spacetime is then given by 



udt 2 + 2dt{ydz + pdx + qdy) — (adx 2 + bdy 2 ), 

(u - 2Vd t K)dt 2 + ZVdtdK + 2dt{pdx + qdy) - {adx 2 + bdy 2 ), 

udt 2 + 2Vdtdr + 2dt(pdx + qdy) - (adx 2 + bdy 2 ), 
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Thus, in the new coordinate system (t,r,x,y), the metric becomes 

/ u V Pr 2 + d x Vr -Qr 2 + d y Vr \ 



(2.20) 
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Similarly, if we take p = VK as new coordinate to replace z, then in the new 
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coordinate system (t,p,x,y), the metric becomes 
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The above procedure can provide a family of solutions depending on 4 coordinates, 
which is relatively complicated. However, for some special cases including more sym- 
metries, the procedure can be much simpler. For instance, if we consider the static 
case with the following ansatz 
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0, u = u(x,y), v = v(x,y), 
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then the solutions can be explicitly obtained. In fact, for this case, G tx , G ty , G zz , G z 
and G zy vanish automatically. By G xx + G yy = 0, we get 



(d 2 x + d 2 v )v = 0, 



(2.25) 



that is to say, v(x,y) is a 2-dimensional harmonic function, which can be easily gener- 
ated by means of complex analytic function f(x + yi). 

Moreover, by G xx — G yy = G xy = G u = and (2.25), we get 



{d 2 x + d 2 y )a = -(\d x a\ 2 + \d y a\ 2 ) + — (\d x v\ 2 + \d y 
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The solution of (2.26) consistent with G xx = G xy = G yy = is given by 

a=^(|^t;| 2 + |a^| 2 ), 



(2.26) 
(2.27) 

(2.28) 



where A is a constant. 

(2.27) is a linear equation of u. Let u = vw, we get 

{d 2 + d 2 )w = —-(d x wd x v + dywdyv). 



(2.29) 



For any given v, the solution of (2.29) can be easily determined by an appropriate 
boundary condition. Here we provide the analytic solution related to v, which is given 
by 

u = v[m + n\nv + (k + j In v)v], (2.30) 
where j, k, m, n are constants, v is the conjugate harmonic function of v, namely, 



v = 9t(f(x + yi)), v = 5s(f(x + yi)), 



(2.31) 
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or 



v = $S(f(x + yi)) 1 v = fft(f(x + yi)). (2.32) 

For other ansatz of the metric, the recursive resolution procedure is basically along 
the same line, that is 

G zz — > (G zx , G zy ) — > (G xy , G xx , Gyy, G tz ) — » ( G t x ) G ty ) — > G tt - (2.33) 



3 Examples 

Example 1. For the case (2.24)~(2.30), if we take the complex function as 

/ = f-§, (3.D 

where Z = x + yi and R is a real number, we have 

sQe 2 + y 2 - R 2 ) y(x 2 + y 2 + fl 2 ) 



Substituting (3.2) into (2.28) and (2.30), we get 

A(x 4 + 2x 2 y 2 + 2R 2 x 2 + y A + i? 4 - 



RJx(x 2 + y 2 - R 2 )(x 2 + y 2 fR 



(3.3) 



This is a more complicated version than that presented in [9]. 

Example 2. For the case (2.24)~(2.30), if we take the complex function as 

/ = cosh ( X -±^-) , (3.5) 



R 



where R is a real number, then we have 



v = dt( f) = cosh — cos — , v = Sf( f) = sinh — sin — . (3.6) 
w; R R K J R R v 7 

Substituting (3.6) into (2.28) and (2.30), we get 



a = ^( cosh2 I- cos2 I)v( cosh l co 4) (3 - 7) 

x y / / x y \ \ 

u = cosh — cos — I m + n In v + (k + j In v) I sinh — sin — j j . (3.8) 
R R V V R RJ J 

Obviously, we can construct infinite solutions in this family by choosing different ana- 
lytic function /. 
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Example 3. For the static case of metric (2.22), if we take the following ansatz 



L = l, iV = sin 2 x, f/ = 0, 



then S — 1 in (2.13) and we can get 



x 



W = 0, P = msinx, Q = /csin x, V = n sin x tarn — , 



(3.9) 



(3.10) 



where fc, m, n, (3 are constants. 

If we take k = a,m = (3 = 0,n = 1 in (3.10), under a suitable coordinate transfor- 
mation, from (3.9) and (3.10) we can get a clear and normal form of metric (2.22) as 
follows 

/ a sin x r cos x \ 
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where a, i? are constants. In this case, the corresponding spacetime has a torus 
structure. 

Example 4. For the dynamic case of metric (2.22), if we take the following 

ansatz 



L = l, iV = sin 2 x, P = Q = d y V=0, 



(3.12) 



then we can get 



or 
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where k, m and n/3 > are constants, and f(t) is an arbitrary function of t. 
In particular, setting k — 1, j3 — |, n = 2 and 



/ = 1 + e sm(cot), 



using (3.15) we get 



V — (1 + £ sin(otf)) sinx * /tan 
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then metric (2.23) becomes 
= diag 
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1 

where e is a constant, (t,p,x,y) are the coordinates. This metric provides a rigorous 
model for the study of the behavior of gravitational waves. 

Example 5. For the periodic case of metric (2.22), if we take the following ansatz 

L=l, N = sm 2 x, V=l, (3.19) 

then we can get 

U = 1, W = 0, a = -2m, (3.20) 

P — fcos(y — h), Q = f cos x sin x sm(y — h), (3.21) 

where m is a constant, and f(t), h(t) are arbitrary functions of t. This metric provides 
a strange rotary gravity. 

Setting / = and making a coordinate transformation 

t = t — r — 2m ln(r — 2m) , x = 8, y = </?, 

metric (2.22) becomes the normal Schwarzschild metric in system (T,r,6,ip): 



[ 2m 
^ = diag j 1 - — , - ( 1 



2m 



2 2 • 2 /) I 

, r , r an v > . 



(3.22) 



Example 6. For metric (2.22), if we take L — N — 1, then we have the following 
solutions with axial symmetry 
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where k, m, (3 are constants, p = y/x 2 + y 2 , f(t) is an arbitrary function, and 
(t, r, x, y) are the coordinates. 

For cases without axial symmetry, we can get solutions as follows: 



/ 4fc 2_ r 2 n 2_ 2 
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k cosh (2x) 
2k sinh (2x) r 
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where k, m, n are constants. 



k cosh (2x) 2k sinh (2x) r —nr 2 ^ 
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4 Conclusion 

All the examples mentioned above always contain some black holes, i.e. singularities 
in Riemann tensor so they are not trivial solutions. These solutions show that 

the canonical metric is really convenient and powerful for solving the vacuum Einstein 
equation, and the procedure of resolution is quite straightforward. 

There are many exact solutions collected in the books like [10, 11]. Most of them 
can be transformed into the canonical form (1.8). Some authors have realized the 
convenience of this kind of form related to the light cone, and the Robinson- Trautman 
type solutions are simple cases of (1.8). However these authors seem to pay attention 
mainly on using the advanced mathematical tools to certain special cases. 

The main characters of our procedure is that: 

(I) It presents a more general framework for getting exact solutions to the Einstein 
equation, which covers most of the known results for empty space. 

(II) The procedure of resolution shows that the Einstein equation of (1.8) is an 
underdetermined system, so that we should take the components a, b as given functions 
determined by coordinate conditions and boundary or initial conditions. 

(III) When we take a, b as given functions, the other components are decoupling 
variables in the system of equations in some sense, so we can use the recursive method 
to solve the Einstein equation. 

(IV) The ansatz (2.3) or (2.24) is natural and necessary for getting exact solutions. 
As to the numerical solution, instead of this ansatz, we should input the appropriate 
values of a, b as boundary or initial conditions to determine other components. 
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